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Effects of QCD Vacuum and Stability of H Dihyperon * 

A. E. Dorokhov, N. I. Kochelev 1 ' 
Joint Institute for Nuclear Research, Laboratory of Theoretical Physics, 140980, Dubna, Moscow Region, Russia 

Within the composite quark model taking into account the interaction of quarks in the bag with 
vacuum fields of QCD the masses of H (I = 0, S = -2, J = 0, l c ) and H* (I = 0, S = -2, J = 1, l c ) 
dihyperons are estimated: Mh = 2.16 GeV, Mh* = 2.34 GeV. It is shown that the leading effect 
giving a stable with respect to strong decays H dihyperon is the instanton interaction forming 
diquarks: q 2 (J = 0, C = 3), q 2 (J = 0, C = 6). With the approach developed the contribution of 
QCD hyperfine interaction is suppressed, and instanton induced three-particle forces in multiquark 
hadrons (q 3 (J = 1/2, C = 8) - channel) are rather small. 

o 

1. Within the quark bag model 0- Q the consideration of multiquark states 0, 0, El was one of the most 
interesting applications. In these calculations the MIT version of bag model was used. As proposed in a static 
' approximation p| , the energy of a multiquark system is determined 0, El] by 

£(*) = £ ! f i + f s * 3 -§+ A ^ a) 

(N ' 



> 



M 2 (R) = E* - <P 2 ) , (^)-E^(f) . ( 2 ) 
(~nJ I flav 
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■ where R is the bag radius, n, is the number of quarks of an i-th type with energy LOi/R (u> Uy d = 2.043 in the ground 
t-H | s-state; m u — nid — 0, m s = 280 MeV), B is an external pressure (B 1 / 4 = 145 MeV), —Zq/R is a contribution of 
zero-mode fluctuations (Zq = 1.84), AE^represents a color-magnetic interaction (a s = 2.2). The stability condition 



O ' dR 



^jjt- = fixes the bag radius. 



In the MIT version of the model the hadron spectrum is specified by the QCD interaction 

Q-T N 
^ ■ 

where N is the total number of quarks, a i (A°) are the spin (color) operators of an z-th quark; puj = /i (rriiR, rrijR) 
. \ determine the strength of a color-spin interaction. In a massless (rrii = 0) case the average of (0 over a hadron state 
^ ■ g™ is expressed through the Casimir operators of spin SU^- | J( J+ 1); color SU^: C3 = 0; spin-color SUq J - defined 
?-h I by quantum numbers of the considered state 
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where J is the total moment. Based on this formula it is possible to formulate the rules |4| analogous to the Hund 
rules from atomic physics. 

For the case of dibaryons rules are the following 5] : The lightest states are those in which quarks are in the most 
symmetrical (antisymmetrical) with respect to color-spin (flavor) representations. 

Such general considerations and calculations based on the MIT model lead to the conclusion [fjj that a flavor-singlet 
six-quark dihyperon H with strangeness —2 and J p = + may be stable with respect to a strong decay. 

Recently, the problem of stability of the H dihyperon has again been discussed after unusual signals from the 
Cygnus X-3 have been registered 9] . As a possible explanation of this effect it was proposed 0] that Cygnus X-3 is a 
star containing the strange matter and it emits H dibaryons with a lifetime th ^ 10 years. The value of H mass is an 
essential point in the determination of its lifetime and confirmation of this hypothesis Tl]. While today experimental 



* Misprints are corrected, literature is updated. Reprinted with partial support from RFBR grant no. 04-02-16445. 
t High Energy Physics Institute, AN Kazakh. SSR 



2 



situation on the emission from Cygnus X-3 stays indefinite the search for multiquark states in cosmic rays and 
on accelerators intensively continues |l3j . That is why correct calculations of multiquark masses, their lifetimes and 
decay modes are so important. In this respect the H dihyperon is most intriguing; object of the investigations. 

The mass of H was estimated within different variants of the bag model pi l?j Il4|. the lattice approach, the QCD 
sum rules method and Skyrme model |15| . The aim of the present work is to calculate the mass of H in the quark 
bag model taking into account the structure of QCD vacuum 0. In [3 such a model was proved to be consistent 
with the method of QCD sum rules and capable of describing the spectroscopy of the ground states of hadrons. 

2. Let us formulate the basic assumptions of the model |l6(. It is known [2| that in the MIT model, it is 
assumed conserning the vacuum structure that in the presence of valence quarks nonperturbative vacuum fully goes 
out of the bag. However, this hypothesis is not compatible with the picture produced by the QCD sum rules 

and it is not self-consistent. Really, the bag constant B characterizing the degree of destroying of vacuum 
-Bmit ~ (130 — 150 MeV) 4 0] is much smaller than the "depth" of nonperturbative vacuum known from the QCD 
sum rules: Eq = — ^ (0 |^G 2 | 0) ~ — (240 MeV) 4 . This means that the vacuum fields practically do not change 
inside the bag-hadron. Therefore, the neglect of the effects of QCD vacuum in the bag becomes physically groundless. 

In works [l6| basic principles were formulated allowing one to consider the QCD vacuum in the bag. 

As a starting point, we take the QCD theory and available information on the behavior of its solution. First, we 
think that in a zero approximation the structure of the solution is such that the low- and high-frequency components 
of the solution are independent of each other. 

Second, we suppose that the (valence) components of the fields with characteristic frequencies to ~ uj q are described 

by solutions of the static bag-model equations (x) — q bag (x) ; At (x) = A T &g (x)j . 

Third, low-frequency (condensate) field components (lu << LJ q ) are assumed to be solutions of the QCD equations 
characterized by a set of numbers: different vacuum condensate quantities ((QQ), (G 2 ),...). Under these assumptions 
the Hamiltonian of the interaction of valence components (q (x) , At (x)) with condensate ones (Q (x) , A vac (x)) is 
restored uniquely through the field transformation: 

* (x) = q (x) +Q(x); A (x) = A T (x) + A vac (x) . (5) 

In additiom there are QCD vacuum fluctuations with w vac >> u> q which may be approximated by the 't Hooft 
interaction [l!j induced by instantons. 

3. In the model ^(J the energy of a hadron is defined as 

M 2 =E 2_ { p2 )i ^- = 0, (6) 

where (P 2 ) — X^Hav n i i K i/R) 2 1S due to the cm. motion of quarks Q and 

E (R) = E kin + AE g + AE vac + AE inst (7) 

is the bag energy. In Ql the kinetic energy of q uarks E'kin and the one-gluon interaction energy AE g are calculated 
as usually in the bag perturbation theory 0, Il8t |2(j : 

2adn = $>ig, (8) 
flav 



A -Eg = °- 11 J Q;s [M 00 + (1 - QA3m s R) M 0s + (1 - 0.25m s R) M ss ] , (9) 

where denotes matrix elements of the operator with respect to spin-color spin states of hadrons. 

As was shown in fTr| a leading contribution to the hadron energy caused by the valence- and condensate-fields 
interaction is generated by the Hamiltonian 

#vac = ^(Q7 <Z + 97°Q)- (10) 
Then by using stationary perturbation theory 

($ \Ht\ *) 

AEvac = ^ ^ I*) =J7 (-OO,0)|$) , (11) 
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U (-oo,0) = ^-T / / dt n T [H vac (h) ...H vac (t n )} 

where | $) is a nonperturbed hadron wave function of the bag model, we have [l6j | 



A^ ac = -„„ ffl^rf _ n ^(QNQ) fl2 (y + °) a y + (12) 

°24 x -l 12 2y(y-l) + a K ' 

tt 2 (0|Hu|0) 2 n5 /~ (x + y)(y + a ) 2 K-l) (OHO) ~ 
+ T^2^o(-o-l) 2 \ M00+ x2[2y (»-!)+ a] (0|l«|0) Ms » + 

%(y + a) 4 x K-l) 2 / (Q|^|Q) \ 2 ~ 1 

>4[2z/(y-i) + «] 2 l(OHO>i "J ' "" 

where (0 0) are quark condensates, y = w s i?, a = m s i?. 

Expressions (|12[1 are absolutely different from BR 3 arising ed hoc in the MIT version. In contrast, in the model 
considered stability of the bag is achieved in a self-consistent manner due to the interaction of quarks with a physical 
vacuum. Moreover, the potential (|12[1 is drastically dependent on the number of quarks with a given flavor, their 

masses and quantum numbers of considered hadron states (the latter is taken into account by the coefficients M [13 )• 
So, the long- wave vacuum fluctuations define the effective quark mass H12(l . At the same time the interaction of 
quarks with a short-wave part of vacuum fluctuations allows us, to a great extent, to explain the mass splitting 
between the terms of SUf(3) hadron multiplets Ha|. Within the model of QCD vacuum as an instanton liquid 



we get the two-particle contribution to the energy [16 

(13) 



= ~ n ^ Vablab jl + (1 + 3~CT a ~a 6 ) j 
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and the contribution of three-particle forces for multiquark states 

p(3) 



A£ W = +noVudsIuds jl + J^ A & (1 + 3~a a ~a b ) + (14) 
- ^d a ^\«\P\] [1 - 3 {^ a ^ b + permutations)] + 

Here n = (0 \g 2 G a ^G^ v \ 0) / (64tt 2 ) is a density of instantons in the model (n = (irp c (~QQ)) 2 /3j, 

Vii...i„ = Q 71 " 2 ^) / ( m nPc--m* ri p c ) , 

2vr 2 

m* = rrii + m*, m* = — (0 |itu| 0) p 2 , 

p c is the characteristic size of an instanton in the QCD vacuum, 

n 

bag i=1 

It should be emphasized that the interaction through instantons (|13|l . (|14|l takes place in a system |) of quarks in 
a zero mode [H El US [H 13 

n 

E(^ + ^)l)=0- (15) 

i=l 
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Diquarks: q 2 (s F , J = 0, 3°) , q 2 h F , J = 1, 6 C ) , triquarks: q 3 (8 F , J = §, 8 C ) , g 3 ( J = | , 10 c ) , etc. are such sys- 
tems. 

We also note that the instanton interaction QlHjl. (|14|1 is consistent only in the first order of perturbation theory 
analogously to that as it was done in the case of an external pion field [25j . 

4. The model parameters. 

As is seen from (Q, the MIT model uses four parameters (B, a s , Zo, m s ). But, the parameter Zq is not well 
grounded, the value of m s is too large, and B poorly agrees with the parameters extracted from the QCD sum rules 
and current algebra. The value a s = 2.2 does not agree with a perturbative expansion in this parameter, which was 
confirmed in one- loop calculations (2(|. In multiquark systems, perturbative calculations with large a s get still more 
uncertain |(J. 

Within the description of energies of the ground states of hadrons it has been proved that it suffices to choose 
the following values of the parameters 

a s = 0.7, m s = 220 MeV, (16) 

and p c = 2 GeV" 1 , (0 \Q t Qi\0) = - (250 MeV) 3 (i = u,d,s) adopted from the model of vacuum 21] and QCD 
sum rules [l7j, respectively. 

5. To calculate the matrix elements of two- and three-particle operators included into AE g , AE vac , AE^L, AE^ t 
it is necessary to know the cluster expansion (dissociation) of the six-quark wave function of H: q e — > q 3 x q 3 , 

„6 _^ „4 w „2 



q —> q x q . The expansion method and wave functions are given in Appendix. 
By using the wave functions (|A.3(I we have for the matrix elements AE g , AE, 



(2) . 



AE? = (-5 M oo ~ 22 MOs + 3 Mss ) /R, (17) 



?H > 1 /H . 67 



AEg = - I —p o - + g-A*~ ) /R, 

2T 

AE^l? = — ^-«o {IudVud + 2I us r) us ) , (18) 

A£ -inst = ^ n {tudVud + 2I us Tj us ) . 

In accordance with the selection rule (|15fl t ne three-particle interaction is nonzero only for the component 
q 3 (I = 0, S = -1, J = 1/2, 70 CJ ). So, with the wave function (|A~3|) we have: 

AE IT = +^-n I udsVuds . (19) 

It should be added that the approximation of coefficients p is given in |2("|: coefficients Ijj in |16| (I s ri s ~ 0.65/o?/o) ; 
/ uds = 0.0244/i? 6 . 

6. By using the above-mentioned relations we obtain the estimation of dihyperon masses 

M H = 2.16 GeV, R H = 5.2 GeV -1 , (20) 

Mh* = 2.34 GeV, H H * = 5.3 GeV" 1 . 

So, our results show that the mass of H is less than 2M\ but above the threshold of NA. Dihyperon H* is absolutely 
unstable: M H * > 2M\. 

In accordance with the estimation of the lifetime of H in AT = 1 weak decays established in work [27]], the state 
with the mass Mh = 2.16 GeV is long-lived: tjj ~ 10 -8 sec. 

Note that in the approach developed a basic cause of the stability of H dihyperon is the interaction of valence 
quarks with short-wave fluctuations and physically is due to the same mechanism by which the mass splitting arises 
in hadron multiplets (tt — p, N — A, and so on splittings). 

We also proved the spectroscopic Hund rule for quark systems At the same time its origin is absolutely 

different. The nonperturbative instanton interaction between a pair of quarks produces strong attraction in a sym- 
metric in color-spin representation and is totally absent in antisymmetric states. The instanton interaction takes 
into account the strong interaction at intermediate distances p c and gives use to the formation of quasibound states, 
diquarks [^Jl^- So, the Hund rule is physically due to the existence of diquarks. 

Moreover, in multiquark system there are multiparticle (n > 2) instanton-induced forces (in colorless baryons such 
interactions are absent because of the selection rule Q15|)). At the same time we show that the contribution of 
three-particle forces to the energy of H state is rather small, AE^ 3 ] t ~ +70 MeV. 

Note also that very recently in the experiment carried by the B.A. Shahbasian group the data that confirm the 
existence of H dihyperon have been obtained [2^] . 
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Appendix. The wave function of dihyperons. 

To construct the wave functions of H and H* dihyperons, we make use of the dissociation method developed in 
[30|. The idea of this method was borrowed from work |(|. That method allows us to construct the wave functions 
of multiquark systems g m g™ with respect to arbitrary dissociation (q mi q ni ) x (q m2 q a2 ) (mi + 7712 = m, n\ + n.2 = n). 
To classify the basis states of multiquarks (q m ), the group Uis is chosen as the group in which the direct production 
of isospin SU% , strangeness Uf, spin ST7/ and color SU% groups is embedded. (Singling out the flavor group SUf is 
not effective because of a strong mixing of the SU3 quantum numbers arising in the multiquark sector.) Thus, the 
quark wave function relative to their quantum numbers transforms by a fundamental representation {1} of the group 
Uis- In this case the scheme of the group reduction of the n-quark-system representation {1™} is the following 

U* n -> U 1S -> SUH° x Ui JC -> ( A. 1 ) 

-» (sul x sui c ) x (C7f x sui c ) -» 

-» su£ x (sul x suf) x erf x (su 2 ; x su£) -> 

-> S't/j x C/f x SU 2 ] x Sirf. 

Using Racah's factorization lemma (1949) and the factorization property of transformation coefficients for direct 
product groups, the transformation of a 6-quark wave function to the dissociation basis may be written in the form 

m 

I (l 3 , l 3 ) l 6 (7, X J C , J rf) , 1»* (S, Xf\ Jsn c s )) J0 C ) = (A.2) 

= E|[ 13 I 1 "" ( J '' A o< « , l"' s (5', A' 5 , J^' s )) A' ; 
('/') 

l 3 (l< (J", A' ', «) , l n * (-5", A£, «)) ^V] 



l 3 ( 1 



n o1"'s^ . I 3 ^l"ol"sjj I (l 3 ; I 3 ) l 6 (l™°l ns )^) 

(I'\' ) ; l n « (/"A ')) I l n ") 1"° (/A ) N 



(V s (5'A'g) ; l n * (5"A'i)) | (l n *, l n *) 1" S (5A S ) 
(A JoMoi A J m ) I (A A ) A (JoAio) 



(A' s J^; X'sJ's^s) I (A' s Ag) A s (JsM5) 
(( «) J', (M) J"| («) Jo, (44') Js (JO) 
(W s ) M ', (/«) M"l (/« Mo, («) A*s (Q£)) • 

Here no (^s) is the number of (u, d) (and s) quarks, 7, S and J are the isospin, strangeness and total spin, respectively, 
the representation of SUq C (SUf) . 

The complete expression of the expansion of the H dihyperon wave function (the H* wave function was given in 
ref. HI) is: 

\H) = 0.867 |ei) - 0.499 |e 2 ) , (A.3) 
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l^{gb(i,2i/2i)g^(iViV2i)— } 



„ | ^,21, 2 21 g_ 2 (l 2 , 3 l 2 ; 2 2l) 



{ Qo G 

2 1 
15 71 

r i 

io vl 

2 

45 

4 | {O-i (i. i 2 , 3 1 2 ; 4 21)} 2 - J± {Q-i (i, i 2 , 3 1 2 ; 2 2i) } 2 



WQ,21, 2 o) g s 2 (l 2 , 3 l 2 ; 2 0)-^ 
{g_! (1, l 2 , 1 2; 2 21) g_! (1, l 2 , 3 l 2 ; 2 2l) + ^} + 

{g^ (i, i 2 , 3 1 2 ; 2 o)} 2 - /J {Q s ; (i, i 2 , 3 1 2 ; 4 o) } 2 4 
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{Q-i (0, 2, 1 1 2 ; 2 21) g_i (0, 2, 3 2; 2 2l) + ^} 



where 



|e 2 ) 



/2 

5^/2 



21 ) Q_ 2 (l 2 , 1 2; 2 21 



^ {g_ x (i, i 2 , 1 2; 2 2i) y - {Q s i (i, i 2 , 3 1 2 ; 2 o)} 



l{g s ;(i,i 2 , 3 iV 



'^{g-^l,! 2 , 3 ! 2 ;^!)} 2 
/ I{fl- 1 (0,2/lV2l)} J 
i{g_ 1 (0,2>2;*2l)} a 
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^{g,^!,! 2 , 3 ! 2 ^!)} 2 
i{(^ 1 (0,2 1 1 lV0)} a 
'i{g_x(0,2»2^2l)} a 



3 o.2 * 



g_ 2 (1 2 , 3 1 2 ; 4 21 



g 3 


^n = 


3;/ = 


2' A °" 


21, Jo 


g 3 


(n = 


2; Jo = 


0,A - 


2, J = 


g 3 


n = 


l;A s = 




= IjMs 



Expression (|A.3|) and the analogous one for H* from |3(| are used to calculate matrix elements of the three-particle 
operator contained in AE^l t . The scalar isospin components with the strangeness —1 : Q-i (0, 2, 1 1 2 ; 2 21) and 
g_i (0, 2, 3 2; 2 2l) only give a nonzero contribution. These components are weighted with the probability 25% in the 
total wave function. 

To compute averages of two-particle operators, one may apply expressions of kind l|A.2|) or to construct the disso- 
ciation q 6 — > q 4 x q 2 . The coefficients of this dissociation may be easily found by using the equations for the Casimir 
operators of the Si//, SUf and SU§ groups. As a result, we have for the basis SU% x SU^ x SU§ 



\H(0 F , J , C ) ) = y ^q 4 (2, 0, 2) q 2 (2, 0, 2) + [q 4 (I 2 , 0, 1 2 ) , 2 (l 2 , 0, 1 2 ) + 
+ g 4 (l 2 ,l,2) g 2 (l 2 ,l,2)+g 4 (2,1,1 2 ) 9 2 (2,1, 1 2 )], 



(A.4) 



\H* (21 F , 1 J , O c )) = ^ 4 (2, 0, 2) q 2 (2, 0, 2) + ^Hg 4 ^ ^ q 2 ^ + (A 5) 

+ \f¥o q4 ^ 2 ) 92 ^ x » 2 ) + /I 94 & i; i2 ) 92 & x ' i2 ) • 

If the dissociations (| A. 3|) - (| A. 5|) are expressed as 

Q 6 = ]>> (^K x ( ( 5 3 ) l = I> (Q 4 ), >< (Q 2 ), > (A.6) 

then the matrix elements of the three- (R3) and two-(i?2) particle operators are calculated with the help 

(Q 6 \R 3 \Q 6 ) =20$>?<Q?|i? 3 |Q?>, (A.7) 



(A. 



where the combinatorial factors take into account the multiple character of interaction in a n-particle antisymmetric 
states. 
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